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We have investigated the Josephson inductance Lj of a Superconducting Cooper Pair Transistor 
(SCPT). We traced the inductance using microwave reflection measurements on a tuned resonance 
circuit in which a SCPT was mounted in parallel to a ^ 200 pH strip line inductance. When the 
inverse of the Josephson inductance, determined on the charge-phase bias plane for a SCPT with 
a Josephson to Coulomb energy ratio of Ej/Ec — 1.75, is integrated over the phase, we obtain a 
current-phase relation, which is strongly non-sinusoidal near the charge degeneracy point. 

PACS numbers: 85.25. Cp, 85.35. Gv, 74.78.Na, 74.25.Fy 



Josephson effect is one of the most spectacular phe- 
nomena in quantum coherent matter. In addition to su- 
perconductors, it has been studied in superfluids^'^ and 
Bose-Einstein condensates^'^. The coupling between the 
ordered states across a Josephson tunnel junction can be 
described by the< Josephson energy E{(^) where de- 
notes the phase difference between the order parameter 
fields on the opposite sides of the junction. Equivalently, 
the nature of a Josephson junction can be characterized 
by its current-phase relation (CPR) = where 

$0 = ^ is the fiux quantum. In a regular tunnel barrier 
junction, the CPR is sinusoidal, which leads to a non- 
linear dynamical inductance of Lj^ = = ^/cCos(/? 
where ^ = ^q^/'^'t^ is the flux corresponding to the phase 
cp^. Injunctions with sinusoidal CPR, \Lj\ is 7r-periodic 
and the inductance changes sign at (p = 7r/2 and at 
(/? = 37r/2. As the inverse inductance describes the en- 
ergy exchange rate of the junction with its surrounding 
electrical circuit, we observe that the junction varies pe- 
riodically between a circuit element either receiving or 
releasing energy, in a sinusoidal fashion in this case. In 
different kinds of junctions, for example, in supercon- 
ducting quantum point contacts (QPC) and in diffusive 
SNS junctions, the CPR can be strongly non-sinusoidal^. 
In the outmost case, like in a QPC with near perfect 
transmission (r ^ 1), I{(p) becomes discontinuous'''^ and 
the region of negative Josephson inductance shrinks to 
zero. Similar discontinuous behavior takes place also in 
two symmetric, classical junctions in series^. 

In a mesoscopic double junction device, charge quanti- 
zation on the island between the junctions starts to play 
a role. This leads to phase fluctuations across the indi- 
vidual junctions, even if the total phase bias acts like a 
classical variable and can be fixed. As a consequence, 
the CPR is strongly modified from the classical result. 
The energy bands, calculable from the Schrodinger equa- 
tion with the Coulombic term acting as the kinetic en- 
ergy operator in the phase picture^'^^, are well known 
and these eigenenergies form the cornerstone of charge 
qubits as well as charge-phase qubits^^ i^^i^^ . A lot of en- 
ergy spectroscopy has been performed on these energy 
levels which, however, does not probe the current-phase 
relation of the ground state unless assumptions on the 



symmetry properties of the bands are made. Here we 
report measurements on the Josephson inductance of a 
superconducting Cooper pair transistor and determine its 
CPR. We present, for the first time to our knowledge, re- 
sults on the current-phase relation that demonstrate the 
strongly non-sinusoidal character near the charge degen- 
eracy point. Our results are based on accurate reflection 
measurements of a charge-phase qubit coupled to a mi- 
crowave resonator. 

In the superconducting state, a tunnel junction stores 
energy in the form of Josephson coupling energy accord- 
ing to E = — Ej cos {(p), where the Josephson energy 
Ej is related to the junction critical current Ic through 
Ic = 2eEj/h. According to electromagnetism, the in- 
ductance is defined by the differential of the total energy 



L 



27r 



(1) 



In the regime of small fluctuations, (p <^1^ Eq. ([T]) yields 
for a single classical Josephson junction (J J) a parametric 
inductance Lj = h/{2el^). For arbitrary phase bias, one 
gets Lj^ = cos (/?. 

Although typically applied to single JJs, these concepts 
can be applied to general junction combinations like a 
single Cooper-pair transistor (SCPT) which consists of 
two JJs and one capacitively coupled gate electrode. The 
full description of SCPT's energy storing capabilities is 
governed by its Hamiltonian^^ 



H = Ec{n 



2Ej cos I cos (9 + 2dEj sin | sin (9, 

(2) 

where Ec = e^/{2C) is the charging energy of the is- 
land, n denotes the number of extra electron charges 
on the island, and = CgVg/e is gate- voltage-induced 
charge in units of single electron on the gate with ca- 
pacitance Cg. n is conjugate to 6/2, where 6 is the su- 
perconducting phase on the island. The asymmetry of 
the two Josephson junctions of the SCPT is described 
hjd= {Eji — Ej2)/2Ej. As we are interested in the 
phase response, the electrostatic energy associated with 
the gate capacitance may be neglected here, contrary to 
the case in Ref. [13 where Josephson capacitance was 
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considered. In order to get the Josephson inductance, 
the eigenenergies of Hamiltonian of Eq. ([2]) have to be 
solved firsifi^ii^. 

Quahtative analysis can be performed in the charge 
basis taking only two charge states into account. Then, 
the Hamiltonian can be written in terms of Pauli spin 



matrices ai 
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^ = - - {B^ax + By&y + Bz^z) , 



(3) 



where the magnetic field components are Bz = 4£^c(l — 
Ug), Bx = 2Ej cos{(p/2) and By = 2dEj cos {(p/2). By 
differentiating the energy eigenvalues, one obtains an an- 
alytic form for the inductance of the ground state 
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^)^^ cos (/9 + 3 + cos 2(p 



, (4) 



where 7 denotes the ratio 7 = Ej/Ec. Because of 
the symmetry in the two level approximation one gets 
Lj| = —Lji for the excited state. When 7 <C 1, the 
two level approximation describes the energy bands ac- 
curately but, when 7^1, more charge states are needed. 
Lji calculated from the two level approximation, com- 
pared to the measured data in Fig. [4] below, is found 
to fall short and numerical calculations including more 
levels are needed to make a detailed comparison. Never- 
theless, Eq. dl]) reveals clearly the main features: 1) the 
supercurrent is strongly suppressed with lowering 7, and 
2) the non-harmonic character grows when approaching 
the charge degeneracy point^^. 

In the experiment, the SCPT's reactive response was 
measured in the L-SET (inductively read single-electron 
transistor) configuration^^, where the junctions are em- 
bedded in an electric LC-oscillator circuit (see also Refs. 
[iilL9.20). The L-SET sample was fabricated using a 
standard e-beam lithography process followed by evap- 
oration of aluminum films at two different angles of inci- 
dence. Circuit diagram of the L-SET sample is shown in 
Fig. [H The SCPT is attached in parallel with an induc- 
tor that is realized using an on-chip microstrip loop. The 
length of the loop is 300 /im, which yields Lm = 168 pH 
inductance. The two parallel capacitors of the resonator 
are made by patterning two large plates ( 1. 3mm x 1.3mm) 
on the substrate (see Fig. 2 of Ref. llJ). The bot- 
tom plate of the capacitor is formed by a 150 nm thick 
niobium ground plane layer under the SisN^ dielectric 
layer whose thickness is around 300 nm. The measured 
resonance frequency yields for the capacitors C = 396 
pF, which agrees well with the parallel plate result using 
Cr = 7.5 for SisN4^. Because our capacitors are non-ideal, 
we have to take into account their ac-losses by introduc- 
ing a shunting resistance whose value R = 1.4 kl^ was 
determined from the Q-factor of the resonator. The cou- 
pling capacitor Cc = 6.5 pF, that decouples the oscilla- 
tor from the driving coaxial cable, is a separate surface 
mount capacitor that is bonded with aluminum wires to 



C= 395.9 pF 
L= 167.8 pH 
ft =1.4 kOhm 
Cn=2.1aF 
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FIG. 1: The circuit diagram of our symmetrized L-SET de- 
vice. Here the SCPT's Josephson inductance Lj is in parallel 
with an on-chip microstrip inductor Lm which form together 
the inductor in the LC-oscillator: 1/L=l/Lj + 1/Lm- Cc is 
the coupling capacitor that decouples the resonator from the 
driving signal generator. 



the top pad of one of the capacitors. The gate modu- 
lation was 2e periodic and from the period we deduce 
Cg = 2.1 aF. From the rf-spectroscopy we conclude that 



Ec/kB = 0.8 K, Ej/kB 
parameter d ^ 0.06. 



1.6 K, and the asymmetry 



The Josephson inductance was measured using a mi- 
crowave reflection measurement setup working around 
fm = 870 MHz. Continuous microwave signal was sent 
from one port of a vector network analyzer through a se- 
ries of attenuators and a circulator to the sample. Part 
of the sent signal was reflected back to the transmission 
line due to the impedance mismatch between Zq = 50 ft 
and the impedance of the resonator circuit, Z (see Fig. 
IB): Vout = rVin, where F = (Z - Zo)/{Z + Zq). The 
reflected signal was then taken from the third port of 
the circulator (isolation 18 dB), amplifled with an ar- 
ray of cold and room temperature ampliflers and con- 
nected to the second channel of the network analyzer. 
The network analyzer gave us the magnitude signal and 
the phase argument (See Figs. 6 and 7 in Appendix I.) of 
the reflection coeflicient which were used to calculate the 
effective impedance of the oscillator circuit. Using stan- 
dard circuit analysis with linear elements, the Josephson 
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FIG. 2: Measured inductance of the ground state on the 
phase bias vs. gate charge plane (inductance scale limited 
to ±250 nH). The plot was obtained using Eq. [5] and the 
data scans presented in the supplementary material. 



inductance then can be calculated from the expression 



FIG. 3: Theoretically calculated Lj of the ground state on 
the Lp vs. rig plane {\Lj\ < ±250 nH). The eigenvalues of the 
SCPT were numerically calculated using 5 charge states. 




Lj = - 

' ' (5) 

where shorthand notations have been assigned to 
impedances of the circuit components Zc = {iu;Cc)~^j 
Zl = iujL, and Zr = {1/R -\- iujC)~^ . The Josephson 
inductance, obtained by using the experimentally deter- 
mined Z, is shown in Fig. O Ideally, the zero phase of 
the reflection signal should be determined separately for 
each bias point by plotting r(/) parametrically on the 
Smith chart and looking for the symmetry axis of the re- 
sulting circle (see, e.g. Ref. [22|). We, however, did this 
check only at one bias point and relied that the refer- 
ence phase will not change substantially. The reflection 
magnitude scale was flxed using reflected signal far away 
from the resonance and the known circuit elements. The 
flnite measurement signal averages the energy bands over 
phase direction. The power used in measurement, -129 
dBm, corresponds to phase amplitude of 0.24 rad (p-p) 
and causes a 3 percent error to the inductance. 

The corresponding inductance plot obtained from the- 
ory is shown in Fig. [31 flve discrete charge states were 
used in the determination of the eigenvalues of Eq. ([2]) . 
The general agreement between Figs. [2] and [3] is remark- 
able, except near the charge degeneracy point. At = 1, 
the gap between the ground and excited states becomes 
small, about A/h = 3.8 GHz, and the upper state can 
either be thermally excited at our effective noise tem- 
perature of ~ 40 mK, governed by the thermalization of 
our circulators, or by multiphoton processes due to the 
microwave drive at 872 MHz. In this case, the effective 
inductance is a mixture of the curvatures of the two en- 
ergy levels making the analysis difficult. The mixing can, 
however, be used to determine the energy level occupa- 
tion in systems where SCPT is considered as a qubiii^. 
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FIG. 4: Comparison of the measured Josephson inductance 
and the theoretical curves at two charge bias values 1.25e 
(filled/blue) and 1.51e (open/green)). The red solid line illus- 
trates Lj obtained from the two level approximation. 



A more quantitative comparison between the experi- 
mental data and the theoretical prediction for Lj{(p) is 
shown in Fig. [H The theory models well the exper- 
imental result when Ej/ks = 1.5 K is employed. The 
SCPT behaves as a non-linear inductor whose inductance 
is tunable by phase bias from a few nH to infinity (in ideal 
case) . There exist substantial regions in the charge-phase 
plane where the Josephson inductance takes negative val- 
ues, and the negative- value regime is diminished as the 
charge degeneracy point is approached. 

The current-phase relation can be obtained by inte- 
grating the inverse of Josephson inductance: The relation 



§^ yields /(^) = 1^ L-/dv + 9{V,) where the 
integration constant giVg) must be set so that the zero 
level of I{ip) is correct. The CPRs at gate charge values 
of = and le are displayed in Fig. [5l We observe that 
far away from the charge degeneracy point, ~ 0, the 
CPR is nearly sinusoidal, while close to = 1 the CPR 
becomes strongly non-sinusoidal. The measured behavior 
at = agrees quite well with theoretical curves calcu- 
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lated using Ej/ks = 1.5 K. The measured variation of 
I{(p) with gate charge is found to be 10% stronger than 
what the theory predicts. We beheve this is partiahy due 
to the uncertainty in the parameter values obtained for 
Ej and Ec using energy level spectroscopy and, more- 
over, to the inaccuracy in the base level of l/Lj, which 
leads to an accumulation of error in the integration pro- 
cess. 

Excitation of the SCPT from the ground state Eq to 
the first excited level Ei is another possible cause for the 
discrepancy between measured and calculated CPRs. We 
note that if the contribution of the upper level is included 
in simple terms {i.e. linear population- weighted addition 
of 1/Lj-values of the ground and excited states), then 
the correction would even worsen the agreement between 
theory and experiment. This is because the contribution 
of the level Ei is smaller than that of the ground state, 
leading to reduced 1/L and its integral. Furthermore, 
inclusion of even higher bands yields similar conclusions 
as their values for l/Lj are also small. 

Effective energy bands with stronger curvature can be 
obtained if one assumes that the total energy is of the 
form 7is{^) = po{(p)Eo{(p) -\- pi{(p)Ei{(p) and that the 
populations po{(p) and pi{(p) depend on cp on the time 
scale of l/fm- Then the derivatives of po{^) and pi{(p) 
may make substantial modifications in the inductance 
and, thereby, in the integrated CPR. For a simple de- 
scription, we assume a thermal distribution at effective 
system temperature T that governs the population of the 
upper state in the driven system according to 



(6) 



where A(n^,(/9) = Ei{ng^(p) — Eo{ng^(p). Below we will 
restrict ourselves only to a qualitative analysis by assum- 
ing that T can be regarded as a constant fitting parame- 
ter over the region of interest, i.e. for the data measured 
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Using this model, we get for the modified CPR Is{y^) 
^^^(bias parameter Ug = le dropped) 



dip 



pie >'bT 
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where I^{(p) is the CPR for the ground state. We have 
fitted Eq. ([7|) to our data in Fig. [5] at = 1 e and 
find an improved agreement when taking T = 0.28 K; 
at = TT this corresponds to pi = 0.30^^. This value 
for T is quite high, well above the cross-over tempera- 
ture of le - 2e periodicity in the gate modulation of our 
SCPT. Consequently, the bath temperature ought to be 
modified strongly by LZ transitions caused by the mea- 
surement drive. Far away from charge degeneracy, LZ- 
transitions are not effective and no modification of the 
CPR is needed in our comparison. 

Non-sinusoidal CPRs have been reported in a number 
of Josephson junctions^. However, in superconducting 
Cooper pair transistors, where Coulomb energy plays an 
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FIG. 5: Current-phase relations (CPR) obtained by integra- 
tion of the inverse Josephson inductance at two gate charge 
values: (nearly sinusoidal) and le (strongly non-sinusoidal). 
Experimental curves are denoted by blue dots while the red 
solid curves display the CPRs calculated from Eq. ([2]); the 
orange solid curve depicts the result calculated from Eq. ([7|) 
using T = 0.28 K. 



important role, the results on (derivatives of) CPRs have 
remained few and unpersuasiv o^^i^^i^^ . Closest to our 
work is the Lj measurement of Ref. but their result 
was distinctly affected by a mixing of the ground and ex- 
cited states due to nonequilibrium quasiparticles. Quasi- 
particles were also a problem in the determination of the 
ground state properties of Nb SCPTs in Ref. where le 
periodicity in Vg was observed contrary to our clean 2e- 
periodic modulation. Moreover, Born and coworkers have 
employed Josephson inductance to characterize qubits^, 
but their ratio of Ej/Ec ~ 30... 2 was too high to make 
observations similar to ours. 

Finally, we believe that our method for constructing 
the CPR will be useful also in studies of other exotic 
Josephson junctions where the critical current is small, 
like S-carbon nanotube-S junctions^^. Typically, criti- 
cal currents in small junctions have been measured by 
applying current bias and recording switching current 
distributions^^. Even though it is not easy to accom- 
plish, it has been recently shown that switching cur- 
rents can be employed to measure current-phase relations 
in break junction samples under circumstances where a 
large Josephson junction is patterned in parallel with the 
break junction^^. Because there is no need for an ex- 
tra junction, we think that our scheme is better suited 
for investigations of novel, fragile junctions where extra 
lithographic steps may be problematic. 

In conclusion, using the phase and magnitude of re- 
fiected microwave signals, we have experimentally mea- 
sured the Josephson inductance in a mesoscopic, super- 
conducting double-tunnel-junction device, i.e., we have 
determined the quantity that is dual to the Josephson 



capacitance 
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We have shown by integrating this 



Josephson inductance that the current-phase relation of a 
SCPT is strongly non-sinusoidal near its charge degener- 
acy point and that substantial tuning by the gate voltage 
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towards nearly sinusoidal CPR at ^ can be reached. 
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I. APPENDIX 

Scans containing the raw data employed in the analy- 
sis: 
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FIG. 6: Magnitude of the reflected signal as a function of the 
phase and charge bias on the superconducting Cooper pair 
transistor. The color scale in dB is given on the right; the full 
reflection corresponds to 9.8 dB. 



phase 




FIG. 7: Phase of the reflected signal in degrees as a function 
of the phase and charge bias. A correction of 25 degrees was 
made to the reference point of the reflection phase in order to 
symmetrize the r(/) on the Smith chart. 
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